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Abstract 

We use Planck data to detect the cross-correlation between the thermal Sunyaev-Zeldovich (tSZ) effect and the infrared emission from the galaxies 
that make up the the cosmic infrared background (CIB). We first perform a stacking analysis towards Planck-conftxmeA galaxy clusters. We detect 
infrared emission produced by dusty galaxies inside these clusters and demonstrate that the infrared emission is about 50 % more extended than the 
tSZ effect. Modelling the emission with a Navarro-Frenk-White profile, we find that the radial profile concentration parameter is C 500 = 1-00*15 j*. 
This indicates that infrared galaxies in the outskirts of clusters have higher infrared flux than cluster-core galaxies. We also study the cross¬ 
correlation between tSZ and CIB anisotropies, following three alternative approaches based on power spectrum analyses: (i) using a catalogue 
of confirmed clusters detected in Planck data; (ii) using an all-sky tSZ map built from Planck frequency maps; and (iii) using cross-spectra 
between Planck frequency maps. With the three different methods, we detect the tSZ-CIB cross-power spectmm at significance levels of (i) 6 cr, 
(ii) 3 cr, and (iii) 4 cr. We model the tSZ-CIB cross-correlation signature and compare predictions with the measurements. The amplitude of the 
cross-correlation relative to the fiducial model is Atsz-cm = 1.2 ± 0.3. This result is consistent with predictions for the tSZ-CIB cross-correlation 
assuming the best-fit cosmological model from Planck 2015 results along with the tSZ and CIB scaling relations. 

Key words, galaxies: clusters - infrared: galaxies - large-scale structure of Universe - methods: data analysis 


1. Introduction 

This paper is one of a set associated with the 2015 release of 
data from the Planck^ mission. It reports the first all-sky de- 
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tection of the cross-correlation between the thermal Sunyaev- 
Zeldovich (tSZ) effect (Sunyaev & Zeldovich 1969, 1972) and 
the cosmic infrared background (CIB; Puget et al. 1996; Fixsen 
et al. 1998; Hauser et al. 1998). An increasing number of ob¬ 
servational studies are measuring the tSZ effect and CIB fluctua¬ 
tions at infrared and submillimetre wavelengths, including inves¬ 
tigations of the CIB with the Spitzer Space Telescope (Lagache 
et al. 2007) and the Herschel Space Observatory (Amblard et al. 
2011; Viero et al. 2012, 2015), and observations of the tSZ effect 
with instruments such as the Atacama Pathfinder Experiment 
(Halverson et al. 2009) and Bolocam (Sayers et al. 2011). In 
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addition, a new generation of CMB experiments can measure 
the tSZ effect and CIB at microwave frequencies (Hincks et al. 
2010; Hall et al. 2010; Dunkley et al. 2011; Zwart et al. 2011; 
Reichardt et al. 2012; Planck Collaboration XXI 2014; Planck 
Collaboration XXX 2014). 

The large frequency coverage of Planck, from 30 to 
857 GHz, makes it sensitive to both of these important probes 
of large-scale structure. At intermediate frequencies, from 70 
to 217 GHz, the sky emission is dominated by the cosmic mi¬ 
crowave background (CMB). At these frequencies, it is pos¬ 
sible to detect galaxy clusters that produce a distortion of the 
CMB blackbody emission through the tSZ effect. At the angu¬ 
lar resolution of Planck, this effect is mainly produced by lo¬ 
cal (z < 1) and massive galaxy clusters in dark matter halos 
(above lO''^ Mq), and it has been used for several studies of clus¬ 
ter physics and cosmology (e.g., Planck Collaboration X 2011; 
Planck Collaboration XI 2011; Planck Collaboration Int. Ill 
2013; Planck Collaboration Int. V 2013; Planck Collaboration 
Int. VIII 2013; Planck Collaboration Int. X 2013; Planck 
Collaboration XX 2014; Planck Collaboration XXIV 2015; 
Planck Collaboration XX 2014; Planck Collaboration XXIV 
2015). At frequencies above 353 GHz, the sky emission is 
dominated by thermal emission, both Galactic and extragalac- 
tic (Planck Collaboration XI 2014; Planck Collaboration XXX 
2014). The dominant extragalactic signal is the thermal infrared 
emission from dust heated by UV radiation from young stars. 
According to our current knowledge of star-formation history, 
the CIB emission has a peak in its redshift distribution between 
z - I and z = 2, and is produced by galaxies in dark matter halos 
of 10"-10'^ Mq; this has been confirmed through the measured 
strong correlation between the CIB and CMB lensing (Planck 
Collaboration XVIII 2014). However, due to the different red- 
shift and mass ranges, the CIB and tSZ distributions have little 
overlap at the angular scales probed by Planck, making this cor¬ 
relation hard to detect. 

Nevertheless, determining the strength of this tSZ-CIB cor¬ 
relation is important for several reasons. Certainly we need 
to know the extent to which tSZ estimates might be contami¬ 
nated by CIB fluctuations, but uncertainty in the correlation also 
degrades our ability to estimate power coming from the “ki¬ 
netic” SZ effect (arising from peculiar motions), which promises 
to probe the reionization epoch (e.g., Mesinger et al. 2012; 
Reichardt et al. 2012; Planck Collaboration Int. XXXVII 2015). 
But as well as this, analysis of the tSZ-CIB correlation allows us 
to better understand the spatial distribution and evolution of star 
formation within massive halos. 

The profile of infrared emission from galaxy clusters is ex¬ 
pected to be less concentrated than the profile of the number 
counts of galaxies. Indeed, core galaxies present reduced in¬ 
frared emission owing to quenching, which occurs after they 
make their first passage inside r ^ Rgoo (Muzzin et al. 2014). 
Using SDSS data, Weinmann et al. (2010) computed the ra¬ 
dial profile of passive galaxies for high-mass galaxy clusters 
(M > lO*"* Mq). They found that the fraction of passive galaxies 
is 70-80 % at the centres and 25-35 % in the outskirts of clusters. 
The detection of infrared emission cluster by cluster is difficult at 
millimetre wavelengths, since the emission is faint and confused 
by the fluctuations of the infrared sky (Galactic thermal dust and 
CIB). Statistical detections of infrared emission in galaxy clus¬ 
ters have been made by stacking large samples of known clus¬ 
ters (Montier & Giard 2005; Giard et al. 2008; Roncarelli et al. 
2010) in IRAS data (Wheelock et al. 1993). The stacking ap¬ 
proach has also been shown to be a powerful method for extract¬ 


ing the tSZ signal from microwave data (e.g.. Lieu et al. 2006; 
Diego (fe Partridge 2009). 

Recently, efforts have been made to model the tSZ-CIB cor¬ 
relation (e.g., Zahn et al. 2012; Addison et al. 2012). Using 
a halo model, it is possible to predict the tSZ-CIB cross¬ 
correlation. The halo model approach allows us to consider dis¬ 
tinct astrophysical emission processes that trace the large-scale 
dark matter density fluctuation, but have different dependencies 
on the host halo mass and redshift. In this paper, we use mod¬ 
els of the tSZ-CIB cross-correlation at galaxy cluster scales. We 
note that the tSZ effect does not possess significant substructure 
on the scale of galaxies, so the tSZ-CIB cross-correlation should 
not possess a shot noise term. 

Current experiments have already provided constraints on 
the tSZ-CIB cross-correlation at low frequencies, between 100 
and 250 GHz. The ACT collaboration set an upper limit p < 0.2 
on the tSZ-CIB cross correlation (Dunkley et al. 2013). George 
et al. (2014), using SPT data and assuming a single correlation 
factor, obtained a tSZ-CIB correlation factor of O.IUqq®; a zero 
correlation is disfavoured at a confidence level of 99 %. 

Our objective in this paper is twofold. First, we characterize 
the CIB emission toward tSZ-detected galaxy clusters by con¬ 
straining the profile and redshift dependence of CIB emission 
from galaxy clusters. Then, we set constraints on the overall 
tSZ-CIB cross-correlation power spectrum, and report the first 
all-sky detection of the tSZ-CIB angular cross-power spectra, at 
a significance level of 4 cr. Our models and results on the tSZ- 
CIB cross-correlation have been used in a companion Planck 
paper Planck Collaboration XXII (2015). 

In the first part of Sect. 2, we explain our modelling approach 
for the tSZ effect and CIB emission at the galaxy cluster scale. 
In the second part of Sect. 2, we describe the model for the tSZ, 
CIB, and tSZ-CIB power and cross-power spectra using a halo 
model. Then in Sect. 3 we present the data sets we have used. 
Sections 4 and 5 present our results for the SED, shape, and 
cross-spectrum of the tSZ-CIB correlation. Finally, in Sect. 6 we 
discuss the results and their consistency with previous analyses. 

Throughout this paper, the tSZ effect intensity will be ex¬ 
pressed in units of Compton parameter, and we will use the best- 
fit cosmology from Planck Collaboration XIII (2015, fourth col¬ 
umn of table 3) using “TT,TE,EE-i-lowP” values as the fiducial 
cosmological model, unless otherwise specified. Thus, we adopt 
T/q = 67.27kms^’ Mpc*, erg = 0.831, and = 0.3156. 

2. Modelling 

To model the cross-correlation between tSZ and CIB 
anisotropies we have to relate the mass, Mgoo, and the redshift, 
z, of a given cluster to tSZ flux, Tgoo, and CIB luminosity Lgoo- 
We define Mgoo (and Rgoo) as the total mass (and radius) for 
which the mean over-density is 500 times the critical density of 
the Universe. Considering that the tSZ signal in the Planck data 
has no significant substructure at galaxy scales, we modelled the 
tSZ-CIB cross-correlation at the galaxy cluster scale. This can be 
considered as a large-scale approximation for the CIB emission, 
and at the Planck angular resolution it agrees with the more re¬ 
fined modelling presented in Planck Collaboration XXX (2014). 

2.1. The thermal Sunyaev-Zeldovich effect 

The tSZ effect is a small-amplitude distortion of the CMB black- 
body spectrum caused by inverse-Compton scattering (see, e.g., 
Rephaeli 1995; Birkinshaw 1999; Carlstrom et al. 2002). Its in¬ 
tensity is related to the integral of the pressure along the line of 
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Table 1. Cosmological and scaling-law parameters for our 
fiducial model, for both the T 500 -AT 500 relation (Planck 
Collaboration XX 2014) and the L 500 -M 500 relation (fitted to 
spectra from Planck Collaboration XXX 2014). 


Planck-SZ cosmology 


. 0.29 ±0.02 

0-8 . 0.77 ± 0.02 

Ho . 67.3 ± 1.4 


Planck-CMB cosmology 

. 0.316 ±0.009 

0-8. 0.831 ±0.013 

Hq . 67.27 ± 0.66 


Afsoo-fsoo 

logT,. -0.19 ±0.02 

087. 1.79 ±0.08 

^sz. 0.66 ± 0.50 


M 500 -L 500 

To . 24.4 ±1.9 

ociB . 0.36 ± 0.05 

jSciB . 1.75 ±0.06 

yciB . 1.70 ±0.02 

Scm . 3.2 ±0.2 

fClB. 1-0 ±0.1 


sight via the Compton parameter, which for a given direction on 
the sky is 

y = r ^^UeTed/. (1) 

J rrieC^ 

Here d/ is the distance along the line of sight, ks, crj, trie, and c 
are the usual physical constants, and tie and are the electron 
number density and the temperature, respectively. 

In units of CMB temperature, the contribution of the tSZ ef¬ 
fect to the submillimetre sky intensity for a given observation 
frequency v is given by 


ATcmb , , 

- = ^(%. (2) 

t CMB 

Neglecting relativistic corrections we have g(v) = xcoth(jr/2) - 
4, with X = hvlik^TcMs). The function g(v) is is equal to 0 at 
about 217 GHz, and is negative at lower frequencies and positive 
at higher frequencies. 

We have used the M 500 - T 500 scaling law presented in Planck 
Collaboration XX (2014), 


Dl(z)Ysoo 

— V 

h 

~2-fQ'sz 

(1 -T)M50o1 

10-4 Mpc^ 

1 * 

0.7 


6x 10'4 Mo 


with E{z) = -\/^m(l + z)^ + for a flat universe. The coef¬ 
ficients Ft, asz, and j3sz are taken from Planck Collaboration 
XX (2014), and are given in Table 1. The mean bias, (1 - b), 
between X-ray mass and the true mass is discussed in detail 
in Planck Collaboration XX (2014, appendix A) and references 
therein. We adopt b = 0.3 here, which, given the chosen cosmo¬ 
logical parameters, allows us to reproduce the tSZ results from 
Planck Collaboration XXII (2015) and Planck Collaboration 
XXIV (2015). 


2.2. Cosmic infrared background emission 

The CIB is the diffuse emission from galaxies integrated 
throughout cosmic history (see, e.g., Hauser & Dwek 2001; 


Lagache et al. 2005), and is thus strongly related to the star- 
formation rate history. The CIB intensity, /(v), at frequency v 
can be written as 


/(v) = 


/ 


. dx(z) j(v,z) 

dz -, 

dz (l±z) 


(4) 


with;^f(^) the comoving distance and j(v,z) the emissivity that is 
related to the star-formation density, pspR, through 


jiy, z) = 


Psfr(z)(1 + z)&ef{(y,z)x^(z) 
K 


(5) 


where K is the Kennicutt (1998) constant (SFR/Ljr = 1.7 x 
10 “'°Mq yr"') and 0efF(v,z) the mean spectral energy distribu¬ 
tion (SED) of infrared galaxies at redshift z. 

To model the T 500 - 4 T 500 relation we use a parametric relation 
proposed by Shang et al. (2012) that relates the CIB flux, L 500 , 
to the mass, M 500 , as follows: 


T50o(v) - To 


4 T 500 


fCIB 


1 X lO'^Mo 


T'(z)0[(l±z)y,7’d(z)], 


( 6 ) 


where Lq is a normalization parameter, 7’d(z) = 7’do(l - 1 - z)“™ 
and 0 [v, Td] is the typical SED of a galaxy that contributes to 
the total CIB emission. 


0 [v, Td] 


v^™Bv(Td), if V < Vo, 
y-rciB, if V > Vo, 


with Vo being the solution of dlog[v^™By(7’d)]/dlog(v) = 
-yciB- We assume a redshift dependence of the form 

'l'(z) = (l+z)‘^™. (7) 


We also define S 500 as 


S 500(v) - 


T50o(v) 

471(1 ± z)xKz)' 


( 8 ) 


The coefficients Tdo^ acm, Pcm, Tcib, and dciB from Planck 
Collaboration XXX (2014) are given in Table 1. We fix the value 
of eciB to 1. In Sect. 2.5, this model of the CIB emission will 
be compared with the Planck measurement of the CIB power 
spectra. We stress that this parametrization can only be consid¬ 
ered accurate at scales where galaxy clusters are not (or only 
marginally) extended. This is typically the case at Planck an¬ 
gular resolution for the low-mass and high-redshift dark matter 
halos that dominate the total CIB emission. 


2.3. Angular power spectra 
2.3.1. The halo model 

To model tSZ, CIB, and tSZ-CIB angular power spectra, we con¬ 
sider the halo-model formalism (see, e.g. Cooray & Sheth 2002) 
and the following general expression 

Lc-Lf +Lf , ( 9 ) 

where A and B stand for tSZ effect or CIB emission, is the 

1-halo contribution, and is the 2-halo term that accounts 

for correlation in the spatial distribution of halos over the sky. 

The 1-halo term is computed using the Eourier trans¬ 

form of the projected profiles of signals A and B weighted by the 
mass function and the A and B emission (see, e.g., Komatsu & 
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t'soo [orcmin^] 

Figure 1. Weight of the mass function as a function of Fsoo- 
Black: the ratio between the number of observed clusters and 
the predicted number of clusters from the Planck-SZ best-fit 
cosmology. Red: parametric formula, Eq. (12), for the selection 
function. 


Seljak 2002, for a derivation of the tSZ angular power spectrum): 


'-e 


— An 




AM- 


d^N 

dMdV 




( 10 ) 


where A^N!AMAV is the dark-matter halo mass function from 
Tinker et al. (2008), dV/dzdQ is the comoving volume element, 
and ^^.ih window functions that account for se¬ 

lection effects and total halo signal. For the tSZ effect, we have 
l^tsz ^ l^N(T 5 oo)T 5 ooyf(M 5 oo,z) and Wn(T50o) is a weight, rang¬ 
ing from 0 to 1 , applied to the mass function to account for the 
effective number of clusters used in our analysis; here Tsoo is the 
tSZ flux, and yt is the Fourier transform of the tSZ profile. For 
the CIB emission we have VTqb = Ssmiv)It{Mso(i,z), where It 
the Fourier transform of the infrared profile (from Eq. 18) and 
‘SsooCv) is given in Eq. ( 8 ). 

The results for the radial analysis in Sect. 4.1.3 show that 
the infrared emission profile can be well approximated by an 
NFW profile (Navarro et al. 1997) with a concentration parame¬ 
ter C 500 = 1 0. The galaxy cluster pressure prohle is modelled by 
a generalized NFW profile (GNFW, Nagai et al. 2007) using the 
best-fit values from Arnaud et al. (2010). 

The contribution of the 2-halo term, accounts for 

large-scale fluctuations in the matter power spectrum that induce 
correlations in the dark-matter halo distribution over the sky. It 
can be computed as 


f^AB,2h 


.4. [Az^lf 
J dzdQ U 


AM- 


d^N 

dMdV 




f 


AM- 


d^N 

dMdV 


W^'^%UM,z)\P{k,z) 


( 11 ) 


(see, e.g., Taburet et al. 2011, and references therein), where 
P{k,z) is the matter power spectrum computed using CLASS 
(Lesgourgues 2011) and hiin(M, z) is the time-dependent linear 
bias factor (see Mo &, White 1996; Komatsu & Kitayama 1999; 
Tinker et al. 2010, for an analytical formula). 

As already stated, at the Planck resolution the tSZ emis¬ 
sion does not have substructures at galaxy scales. Consequently 


there is no shot-noise in the tSZ auto-spectrum and all the tSZ- 
related cross-spectra. Considering the method we used to com¬ 
pute the CIB auto-correlation power spectra, the total amplitude 
of the 1 -halo term should include this “shot-noise” (galaxy auto¬ 
correlation). We have verified that, at the resolution of Planck, 
there is no signihcant difference between our modelling and a di¬ 
rect computation of the shot-noise using the sub-halo mass func¬ 
tion from Tinker & Wetzel (2010), by comparing our modelling 
with Planck measurements of the CIB auto-spectra. 

2.3.2. Weighted mass-function for selected tSZ sample 

Some of our analyses are based on a sample selected from a tSZ 
catalogue. However, we only consider confirmed galaxy clus¬ 
ters with known redshifts. Therefore, it is not possible to use 
the selection function of the catalogue, which includes some un¬ 
confirmed clusters. To account for our selection, we introduce a 
weight function, Wn, which we estimate by computing the tSZ 
flux T 500 as a function of the mass and redshift of the clusters 
through the Y-M scaling relation. Then we compute the ratio 
between the number of clusters in our selected sample and the 
predicted number (derived from the mass function) as a function 
of the flux T 500 . We convolve the mass function with the scatter 
of the Y-M scaling relation, to express observed and predicted 
quantities in a comparable form. Uncertainties are obtained as¬ 
suming a Poissonian number count for the clusters in each bin. 
Finally, we approximate this ratio with a parametric formula: 

Wn = erf(660 Tsoo - 0.30). (12) 

This formula is a good approximation for detected galaxy clus¬ 
ters, with isoo > 10 arcmin^. The weight function is presented 
in Fig. 1 . This weight applied to the mass function is degenerate 
with the cosmological parameters, and thus cancels cosmologi¬ 
cal parameter dependencies of the tSZ-CIB cross-correlation. 

Moreover, considering that the detection methods depend on 
both T 500 and 6500 , a more accurate weight function could be de¬ 
fined in the Tsoo-^soo plane and convolved with the variation of 
noise amplitude across the sky. However, given the low number 
of clusters in our sample, we choose to define our weight only 
with respect to Tsoo- 

2.4. Predicted tSZ-CIB angular cross-power spectrum 

In Fig. 2 we present the predicted tSZ-CIB angular cross-power 
spectra from 100 to 857 GHz, for the fiducial cosmological 
model and scaling-relation parameters listed in Table 1. The tSZ 
angular auto-spectram is dominated by the 1 -halo term, while 
the CIB auto-spectrum is dominated by the 2-haIo-term up to 
I ^ 2000. Thus we need to consider both contributions for the 
total cross-power spectrum. 

At low I, we observe that the 2-halo term has a similar ampli¬ 
tude to the 1-halo term at all frequencies. The 1-halo term com¬ 
pletely dominates the total angular cross-power spectrum up to 
I ^ 2000. We also notice that the cross-power spectrum is highly 
sensitive to the parameters ^cib and eciB- Indeed, these two pa¬ 
rameters set the overlap of the tSZ and CIB window functions in 
mass and redshift. Similarly, the relative amplitude of the 1-halo 
and 2 -halo term is directly set by these parameters. 

In Fig. 3, we present the redshift and mass distribution of 
tSZ and CIB power. These distributions are different for differ¬ 
ent multipoles. Given the angular scale probed by Planck, we 
show them for the specihc multipole € - 1000. We notice that 
the correlation between tSZ and CIB, at a given frequency, is 
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Figure 2. Predicted tSZ-CIB cross-correlation from 100 to 857 GHz for the fiducial model, where the tSZ signal is expressed in 
Compton parameter units. The blue dashed line presents the prediction for the 1-halo term, the green dashed line for the 2-halo term 
and the red solid line for the total model. 


determined by the overlap of these distributions. Clusters that 
constitute the main contribution to the total tSZ power are at low 
redshifts (z < 1). The galaxies that produce the CIB are at higher 
redshifts (1 < z < 4). The mass distribution of CIB power peaks 
near M 500 = 10'^ Mq, while the tSZ effect is produced mostly 
by halos in the range 10 '"^ Mq < M 500 < IO'^Mq. 

We dehne the correlation factor between tSZ and CIB sig¬ 
nals, p, as 

(^tSZ-CIB 

^-G77- (13) 

/(^tSZ-tSZ (^CIB -CIB y' ^ 

Figure 4 shows that we derive a correlation factor ranging from 
0.05 to 0.30 at Planck frequencies. This agrees with the values 
reported from other tSZ-CIB modelling in the literature (Zahn 
et al. 2012; Addison et al. 2012), ranging from 0.02 to 0.34 at 95, 
150, and 220 GHz. The difference in redshift and mass distribu¬ 
tions of tSZ and CIB signals explains this relatively low degree 
of correlation. 

We also observe that the tSZ-CIB correlation has a minimum 
around { - 300, and significantly increases at higher multipoles. 
At those multipoles, the tSZ effect is dominated by low-mass and 
higher-redshift objects, overlapping better with the CIB range of 
masses and redshifts, which explains the increase of the corre¬ 
lation factor. The frequency dependence of the tSZ-CIB correla¬ 
tion factor can be explained by the variation of the CIB window 
in redshift as a function of frequency. At high frequencies, we 
observe low-redshift objects (with respect to other frequencies, 
but high-redshift objects from a tSZ perspective). On the other 
hand, at low frequencies, we are sensitive to higher redshift, as 
shown in Fig. 3. 

2.5. Comparison with tSZ and CIB auto-spectra 

We fixed Tdo, ctcib, ySciB, Tcib, and dciB to the values from Planck 
Collaboration XXX (2014). We hx eciB to a value of 1.0, and we 
ht for Lo in the multipole range 100 < ^ < 1000 using CIB 
spectra from 217 to 857 GHz. We notice that the value of ecm is 


closely related to halo occupation distribution power-law index 
and highly degenerate with Qm. 

In Fig. 5, we compare our modelling of the tSZ and CIB 
spectra with measured spectra from the Planck tSZ analysis 
(Planck Collaboration XXII 2015) and Planck data at 217, 
353, 545, and 857 GHz (Planck Collaboration XXX 2014). For 
more details of these measurements see the related Planck pa¬ 
pers referenced in Planck Collaboration I (2014) and Planck 
Collaboration I (2015). Here, we used the Planck-SZ best-fit 
cosmology presented in Table 1. We observe that our model re¬ 
produces the observed auto-power spectra for both tSZ and CIB 
anisotropies, except at low £ (below 100), where the CIB power 
spectra are still contaminated by foreground Galactic emission. 
For this reason, in this £ range the measured CIB power spec¬ 
tra have to be considered as upper limits. The figures show the 
consistency of the present CIB modelling at cluster scale with 
the modelling presented in Planck Collaboration XXX (2014) in 
the multipole range covered by the Planck data. Note that the 
flatness of the 1-halo term for CIB spectra ensures that this term 
encompasses the shot-noise part. 


3. The data 

3.1. Planck frequency maps 

In this analysis, we use the Planck full-mission data set (Planck 
Collaboration I 2015; Planck Collaboration VIII 2015). We 
consider intensity maps at frequencies from 30 to 857 GHz, 
with V.l pixels, to appropriately sample the resolution of the 
higher-frequency maps. For the tSZ transmission in Planck 
spectral bandpasses, we use the values provided in Planck 
Collaboration IX (2014). We also used the bandpasses from 
Planck Collaboration IX (2014) to compute the CIB transmis¬ 
sion in Planck channels for the SED. For power spectra analy¬ 
ses, we use Planck beams from Planck Collaboration IV (2015) 
and Planck Collaboration VII (2014). 
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-^500 [^o] 

Figure 3. Top: predicted distribution of the tSZ and CIB power 
as a function of the redshift at ^ = 1000. Bottom: predicted dis¬ 
tribution of the tSZ and CIB power as a function of the host 
halo mass at ^ = 1000. The black dashed line is for the tSZ ef¬ 
fect, while the dark blue, light blue, green, yellow, orange, and 
red dashed lines are for CIB at 100, 143, 217, 353, 545, and 
857 GHz respectively. The vertical solid black line shows the 
maximum redshift in PSZ2 (top panel) and the minimal M 500 in 
PSZ (bottom panel). 


3.2. The Planck SZ sample 

In order to extract the tSZ-CIB cross-correlation, we search for 
infrared emission in the direction of clusters detected through 
their tSZ signal. In this analysis, we use galaxy clusters from the 
Planck SZ catalogue (Planck Collaboration XXVII 2015, PSZ2 
hereafter) that have measured redshifts. We restrict our analysis 
to the sample of confirmed clusters to avoid contamination by 
false detections (see Aghanim et al. 2014, for more details). This 
leads to a sample of 1093 galaxy clusters with a mean redshift 
z ^ 0.25. 

From this sample of clusters, we have built a reprojected 
tSZ map. We use a pressure profile from Arnaud et al. (2010), 
with the scaling relation presented in Planck Collaboration XX 
(2014), as well as the size ( 6500 ) and flux (T 500 ) computed from 
the 2-D posterior distributions delivered in Planck Collaboration 
XXVII (2015).^ We project each cluster onto an oversampled 
grid with a pixel size of 0.1 x ^500 (using drizzling to avoid flux 
loss during the projection). Then we convolve the oversampled 

^ http://pla.esac.esa.int/pla/ 



Figure 4. Predicted correlation factor of the tSZ-CIB cross¬ 
spectrum from 100 to 857 GHz. The grey shaded area represents 
the range of values predicted for p from Zahn et al. (2012) for 
various models at 95, 150, and 220 GHz. 


map with a beam of 10' FWHM. We reproject the oversampled 
map onto a HEALPix (Gorski et al. 2005) full-sky map with 117 
pixels (Aside = 2048) using nearest-neighbour interpolation. 

3.3. IRAS data 

We use the reprocessed IRAS maps, IRIS (Improved 
Reprocessing of the IRAS Survey, Miville-Deschenes & 
Lagache 2005) in the HEALPix pixelization scheme. These offer 
improved calibration, zero level estimation, zodiacal light sub¬ 
traction, and destriping of the IRAS data. The IRIS 100, 60, and 
25 pm maps are used at their original resolution. Missing pixels 
in the IRIS maps have been set to zero for this analysis. 

4. Results for tSZ-detected galaxy clusters. 

In this section, we present a detection of the tSZ-CIB cross¬ 
correlation using known galaxy clusters detected via the tSZ ef¬ 
fect in the Planck data. In Sect. 4.1, we focus on the study of 
the shape and the SED of the infrared emission towards galaxy 
clusters. Then Sect. 4.2 is dedicated to the study of the tSZ-CIB 
cross-power spectrum for confirmed tSZ clusters. 

4. 1. Infrared emission from clusters 

4.1.1. Stacking of Planck frequency maps 

To increase the significance of the detection of infrared emission 
at galaxy-cluster scales, we perform a stacking analysis of the 
sample of SZ clusters defined in Sect. 3. Following the methods 
presented in Hurier et al. (2014), we extract individual patches 
of 4° x4° from the full-sky Planck intensity maps and IRIS maps 
centred at the position of each cluster. The individual patches are 
re-projected using a nearest-neighbour interpolation on a grid of 
0'.2 pixels in gnomonic projection to conserves the cluster flux. 
We then produce one stacked patch for each frequency. To do so, 
the individual patches per frequency are co-added with a con¬ 
stant weight. This choice accounts for the fact that the main con¬ 
tribution to the noise, i.e., the CMB, is similar from one patch to 
another. Furthermore, it avoids cases where a particular cluster 
dominates the stacked signal. Considering that Galactic thermal 
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Figures. Upper right panel: observed tSZ power spectrum: Planck data from Planck Collaboration XXI (2014) (black symbols), 
ACT data (Reichardt et al. 2012) (light blue symbols), and SPT data (Sievers et al. 2013) (orange symbols); with our hducial model 
(dashed blue line). Other panels: observed CIB power spectra, with Planck data from Planck Collaboration XXX (2014) (black) 
and our hducial model (dashed blue line). These panels show auto- and cross-power spectra at 217, 353, 545, and 857 GHz. The 
dotted blue lines show the 1-halo plus shot-noise term for our hducial model. 


dust emission is not correlated with extragalactic objects, emis¬ 
sion from our Galaxy should not bias our stacking analysis. We 
verihed that the stacking is not sensitive to specihc orientations, 
which may be produced by the thermal dust emission from the 
Galactic plane. We produced a stacked patch per frequency for 
the whole cluster sample, and for two large redshift bins (below 
and above z - 0.15). 


In Fig. 6, we present the stacked signal at the positions of 
the sample of conhrmed SZ clusters in Planck data from 30 to 
857 GHz and in IRIS maps from 100 to 25 pm. At low frequen¬ 
cies (below 217 GHz) we observe the typical tSZ intensity decre¬ 
ment. However, at 353 and 545 GHz we see a mix of the positive 
tSZ signal and infrared emission. We also note that the infrared 
emission can also be observed at 217 GHz where the tSZ effect is 
negligible. We note the presence of signihcant infrared emission 
in the Planck 857 GHz channel, where the tSZ signal is negligi¬ 
ble. Similarly we hnd a signihcant infrared signal in the IRAS 
100 and 6Qpm bands. 


4.1.2. The SED of galaxy clusters 

Each stacked map from 70 GHz upwards is created at a resolu¬ 
tion of FWHM = 13'.^ We measure the Hux, F{v), in the stacked 
patches from 30 to 853 GHz through aperture photometry within 
a radius of 20' and compute the mean signal in annuli ranging 
from 30' to 60' to estimate the surrounding background level. We 
use aperture photometry in order to obtain a model-independent 
estimation of the total hux, without assuming a particular shape 
for the galaxy cluster prohle. Thus, the hux is computed as 


F{v) = K, 




N, 


/•< 20 ' 


V/-<20' 


30’<r<60’ 


-^30'<r<60' 


AQ, (14) 


where Ayj, is the pixel p separated from the centre of the map by 
a distance r in the stacked map Ay, AQ is the solid angle of one 
pixel, Nx is the number of pixels that follow the condition X, and 
Ky is the bias for the aperture photometry (equal to one except at 
30 and 44 GHz). The 30 and 44 GHz channels have a lower an¬ 
gular resolution and thus the aperture photometry does not mea¬ 
sure all the signal. We compute the factor Ky by assuming that 

^ At 30 GHz and 44 GHz, we keep the native resolution of the inten¬ 
sity maps, i.e., FWHM values of 32134 and 27112, respectively. 
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Figure 6. From left to right and top to bottom: observed stacked 
intensity maps at 30, 44, 70, 100, 143, 217, 353, 545, 857, 3000, 
5000, and 12 000 GHz at the positions of confirmed SZ clus¬ 
ters. All maps are at the native angular resolution of the rele¬ 
vant Planck channel. The measured intensity increases from dark 
blue to red. 


the physical signal at 143 GHz has the same spatial distribution 
as that at 30 or 44 GHz: 


Ky 


F(143) 

F*(143,v)’ 


(15) 


where the flux ^*(143, v), is computed on the 143 GHz stacked 
map after being set to the resolution of the 30 and 44 GHz chan¬ 
nels. 

From the scaling laws presented in Sect. 2, it is possible to 
predict the expected SED of a cluster. The total tSZ plus infrared 
emission measured through stacking of Ad clusters can be writ¬ 
ten as 

A'd 

F(y)^J][8(y)Y‘^ + S‘,,o(y)], (16) 

i 

where S^oaiy) - a L5oo(y)l4nx^(z) is based on values from 
Table 1 and we fit for Lq (see Eq. 6 for the L 500 expression), 
which sets the global amplitude of the infrared emission in clus¬ 
ters. 

Eigure 7 presents the derived SED towards galaxy clusters 
compared to the tSZ-only SED. The observed flux at high fre¬ 
quencies, from 353 to 857 GHz, calls for an extra infrared com¬ 
ponent to account for the observed emission. We also present in 
Eig. 8 the same SED for two wide redshift bins, below and above 
z - 0.15, with median redshift 0.12 and 0.34, respectively. We 
observe that most of the infrared emission is produced by objects 



0 200 400 600 800 


u [GHz] 

Figure 7. Observed SED of the stacked signal towards galaxy 
clusters, from 30 to 857 GHz. In blue we show the tSZ contri¬ 
bution to the total SED and in red the total SED considering 
infrared emission. 


at z > 0.15. We compare this stacking analysis to the SED pre¬ 
diction from the scaling relation used to reproduce results pre¬ 
sented in Eig. 7 (red lines). This shows that the modelling repro¬ 
duces the observed redshift dependence of the infrared flux of 
galaxy clusters. 

Uncertainties on F{y) are mainly produced by contamination 
from other astrophysical components. To estimate the induced 
contamination level, we extract fluxes, F'at 1000 random po¬ 
sitions, q, across the sky with the same aperture photometry as 
the one presented in Sect. 4.1.2. To derive a realistic estimation 
of the noise, we avoid the Galactic plane area for the random 
positions, since this area is not represented in our cluster sample 
(for details of the sky coverage of the PSZl and PSZ2 catalogues 
see Planck Collaboration XXIX 2014 and Planck Collaboration 
XXVII 2015). In the stacking process each cluster is consid¬ 
ered uncorrelated with the others. Indeed, considering the small 
number of objects, we can safely neglect correlations induced by 
clustering in the galaxy cluster distribution. The Fy uncertainty 
correlation matrix is presented in Pig. 9. It only accounts for un¬ 
certainties produced by uncorrelated components (with respect 
to the tSZ effect) in the flux estimation. Prom 44 to 217 GHz, 
the CMB anisotropies are the main source of uncertainties. This 
explains the high level of correlation in the estimated fluxes. The 
30 GHz channel has a lower level of correlation due to its high 
noise level. At higher frequencies, from 353 to 857 GHz, dust 
residuals becomes the dominant source in the total uncertainty, 
which explains the low level of correlation with low frequency 
channels. Contamination by radio sources inside galaxy clusters 
can be neglected (at least statistically), since the measured fluxes 
at 30 and 44 GHz agree with the tSZ SED. 

4.1.3. Comparison between tSZ, IR, and galaxy-number 
radial profiles 

We computed the radial profile of each stack map from 70 to 857 
GHz in native angular resolution.These profiles were calculated 
on a regular radial grid of annuli with bins of width Ar - 1', 
allowing us to sample the stacked map at a resolution similar 
to the Planck pixel size. The profile value in a bin is defined as 
the mean of the values of all pixels falling in each annulus. We 
subtract a background offset from the maps prior to the profile 
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Figures. Observed SED of the stacked signal toward galaxy 
clusters, from 30 to 857 GHz. Top: objects below z = 0.15. 
Bottom: objects above z = 0.15. In blue we show the tSZ con¬ 
tribution to the total SED and in red the total SED considering 
infrared emission. 


computation. The offset value is estimated from the surrounding 
region of each cluster (30' < r < 60'). The uncertainty associ¬ 
ated with this baseline offset subtraction is propagated into the 
uncertainty of each bin of the radial profile. 



R [orcmin] 



R' [orcmin] 



Figure 9. Correlation matrix of F{v) estimated from 1000 ran¬ 
dom positions across the sky. 


Figure 10. Top: observed radial profile of the stacked signal to¬ 
ward galaxy clusters at the native angular resolution at 70 GHz 
(in dark blue), 100 GHz (blue), 143 GHz (light blue), 217 GHz 
(green), 353 GHz (orange), 545 GHz (red), 857 GHz (dark red), 
and 100 pm (black). Bottom: same as top panel, but as a function 
of the rescaled radius R' - RifivISm). 


In Fig. 10, we present the profile normalized to one at the 
centre. We observe that profiles derived from low-frequency 
maps show a larger extension due to beam dilution. The smallest 
extension of the signal is obtained for 353 GHz, then it increases 
with frequency. For comparison we also display the profile at 
100jum, which shows the same extension as the profiles at 217 
and 857 GHz where there is no significant tSZ emission. For il¬ 
lustration, in the bottom panel of Fig. 10, we display the profiles 
as a function of the rescaled radius R' - RiSvISgsi), with Sv 
the FWHM of the beam at frequency v. Under the assumption 
that the tSZ profile is Gaussian, this figure allows us to directly 
compare the extension of the signal at all frequencies; it illus¬ 
trates the increase of the signal extension with frequency, except 
for the 217 GHz profiles, which have the same size as the high- 
frequency profiles. 
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We then define the extension of the Planck profiles as £(v), 
where 


&Hy) 


= 471 In 2 


7 rp(r,v)drf 
. f P(r, v)dr , 


-St 


(17) 


with p(r, v) the profile of the stacked signal at frequency v. 
Integration of the profiles is performed up to r = 30'. The quan¬ 
tity £(v) is equivalent to a FWHM for a Gaussian profile; how¬ 
ever, we notice that the profile of the stacked signal deviates from 
a Gaussian at large radii. This increases the values obtained for 
the profile extension. We estimate the uncertainty on the profile 
extension using 1000 random positions on the sky. 



V [GHz] 


Figure 11. Variation of the spatial extension, £(v), of the stacked 
signal from 70 to 3000 GHz (100/rm). The blue line shows the 
expected value for the GNFW profile from Arnaud et al. (2010) 
and 0500 values from PSZl, the orange line shows the value we 
derive with the Xia et al. (2012) profile, and the red line shows 
the value for an NFW profile with C 500 = 1 0. 


In Fig. 11, we present the variation of £(v) as a function of 
frequency from 70 GHz to 3000 GHz (100/tm). We see lower 
values of £(v) at low frequencies. The observed signal is com¬ 
posed of two separate components, the tSZ effect and the in¬ 
frared emission from clusters. The variation of &{v) is produced 
by the difference between the spatial extension of the tSZ effect 
and the infrared emission. At frequencies dominated by the tSZ 
signal (from 70 to 143 GHz), we observe £(v) = 13(6 + 0(1. 
The expected value for the Arnaud et al. (2010) GNFW profile 
and 0500 values for the Planck clusters is £(v) = 13(6. We ob¬ 
serve £(v) = 20(0 + 0(5 at 217, 857, and 3000 GHz. At 217 GHz, 
where the tSZ signal is almost null, £(v) is dominated by the 
infrared emission and is similar to the signal found at high fre¬ 
quencies (857 GHz and 100jum). Considering an NFW profile 
for infrared emission, 

p{r) cx -—i--, (18) 

V500^)(l +C500^) 

and 0500 values for the Planck clusters, the previous result trans¬ 
lates into constraints on C500, giving C500 = l-OO);^}^. For com¬ 
parison, we also display in Fig. 11 the prediction based on the 
profile used in Xia et al. (2012), which assumes a concentration 


9 lM\ 


where M* is the mass for which v(M, z) — 6c I {p^criM'^ is equal 
to 1 (Bullock et al. 2001), with 5c the critical over-density, Dg the 
linear growth factor, and cr{M) the present-day rms mass fluctu¬ 
ation. 

This model leads to £(v) = 17(6. As a consequence, our re¬ 
sults demonstrate that galaxies in the outskirts of clusters give 
a larger contribution to the total infrared flux than galaxies in 
the cluster cores. Indeed, star formation in outlying galaxies 
is not yet completely quenched, while galaxies in the core no 
longer have a significant star formation rate. This result is con¬ 
sistent with previous analyses that found radial dependence for 
star-forming galaxies (e.g., Weinmann et al. 2010; Braglia et al. 
201 1; Coppin et al. 2011; Santos et al. 2013; Muzzin et al. 2014). 

We also model the radial dependence of the average specific 
star-formation rate (SSFR) as |^1 - Aqexp 7aqr//?5oo)], where 
1 - Aq is the ratio between the SSFR of a core galaxy and 
an outlying galaxy and Oq is the radial dependence of the in¬ 
frared emission suppression. We adopt the profile from Xia et al. 
(2012) for the galaxy distribution and, using Aq ^ 0.7, we derive 
ry = 0 5+°-^ 


4.2. Cross-correlation between the tSZ catalogue and 
temperature maps 

4.2.1. Methodology 


We focus on the detection of the correlation between tSZ and 
infrared emission at the positions of confirmed galaxy clusters. 
To do so, we use a map constructed from the projection of 
confirmed SZ clusters on the sky, hereafter called the “repro¬ 
jected tSZ map,” (see end of Sect. 3.2). We measure the tSZ-CIB 
cross-correlation by computing the angular cross-power spec¬ 
trum, C^”, between the reprojected tSZ map, yc, and the Planck 
intensity maps, Ty, from 100 GHz to 857 GHz. We note that yc 
is only a fraction of the total tSZ emission of the sky, y. The 
cross-spectra are 

= g{v)Cp + Cf (20) 

where c^'TiB(y) cross-correlation between the reprojected 
tSZ map (in Compton parameter units) and the CIB at frequency 
V. Considering that the tSZ power spectrum is dominated by the 
1-halo term, we have Cp ^ We mask the thermal dust 

emission from the Galaxy, keeping only the cleanest 40 % of the 
sky. This mask is computed by thresholding the 857 GHz Planck 
full sky map at 30' FWHM resolution. We verified that we derive 
compatible results with 30 and 50 % of the sky. We bin the cross¬ 
power spectra and correct them for the beam and mask effects. 
In practice, the mixing between multipoles induced by the mask 
is corrected by inverting the mixing matrix M*/, between bins of 
multipoles (see Tristram et al. 2005). 

We estimate the uncertainties on the tSZ-CIB cross-spectra 
as 


II 

(N 

G 

(c7’' 


V ‘ ^ (2^+l)/sky 

V t 



( 21 ) 


We stress that due to cosmic variance the uncertainties on the 
spectra are highly correlated from frequency to frequency. 
We see in Fig. 12 that cross-spectra at different frequencies show 
similar features. The covariance matrix between cross-spectra at 
frequencies v and v' can be expressed as 

j (2^+l)/sky^ ^ ^ ( e J- 


( 22 ) 


-0.13 


(19) 
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Figure 12. Observed cross-correlation between the tSZ map (jc) of known clusters and the Planck frequency maps {Ty), from 100 
to 857 GHz. The data are presented as black circles, the tSZ auto-correlation is in green and the total model y^-Ty accounting for 
tSZ-CIB correlation is in red for the best-fit of the CIB and tSZ spectra. All power spectra are presented in Compton-parameter 
units. Uncertainties are dominated by foreground residuals; thus they are highly correlated from one channel to another. 



Then we propagate the uncertainties through the bins and the 
mixing matrix M^^*. We verify using Monte Carlo simulations 
that we derive compatible levels of uncertainty. In the following 
analysis, we use the full covariance matrix between multipole 
bins. 

We also consider uncertainties produced by the Planck band- 
passes (Planck Collaboration IX 2014). This source of uncer¬ 
tainty reaches up to 20 % for the tSZ transmission at 217 GHz. 
We also account for relative calibration uncertainties (Planck 
Collaboration VIII 2015) ranging from 0.1 % to 5 % for differ¬ 
ent frequencies. We verify that our methodology is not biased 
by systematic effects by cross-correlating Planck intensity maps 
with nominal cluster centres randomly placed on the sky, and we 
observe a cross-correlation signal compatible with zero. 

4.2.2. Results 

In Fig. 12, we present the measured angular cross-power spec¬ 
tra and our fiducial model. For convenience, all spectra are 
displayed in Compton-parameter units; as a consequence the 
tSZ auto-correlation has the same amplitude at all frequencies. 
We note that at 217 GHz the tSZ transmission is very faint. 
Thus, it induces large uncertainties when displayed in Compton- 
parameter units; at this specific frequency, the cross-power spec¬ 
trum is dominated by the tSZ-CIB contribution, since the CIB 
emission dominates over the tSZ emission at 217 GHz (g(v) be¬ 
comes negligible in Eq. 20). Uncertainties are highly correlated 
from one channel to another, explaining the similar features we 
observed in the noise at all frequencies. 

In order to address the significance of the tSZ-CIB corre¬ 
lation in the measured cross-spectra we consider three cases to 
describe the angular cross-power spectra and we compute the;^^ 
for each case and each frequency: 

- in Case 1, no tSZ auto-correlation and no tSZ-CIB correla¬ 
tion; 


- in Case 2, only a tSZ auto-correlation contribution; 

- in Case 3, both tSZ and tSZ-CIB spectra. 

We present the derivedvalues in Table 2. 

At low frequencies, 100 and 143 GHz, the measured signal 
is completely dominated by the tSZ auto-spectrum. The con¬ 
sistency between the observed spectrum and the predicted one 
demonstrates that fluxes from the Planck SZ clusters are con¬ 
sistent with our measurement. At intermediate multipoles the 
tSZ-CIB contribution has a similar amplitude to the contribution 
of tSZ auto-spectra. At 217 and 353 GHz we observe a higher 
value for the in Case 3 compared to the value for in Case 2. 
However, the difference between these values is not significant, 
considering the number of degree of freedom per spectra. Thus 
the tSZ-CIB contribution is not significant for these frequencies. 
However, at 545 and 857 GHz we detect a significant excess with 
respect to the tSZ auto-correlation contribution. We detect the 
tSZ-CIB contribution for low-redshift objects at 5.8 and 6.0 cr at 
545 and 857 GHz, respectively. 

Table 2. Value of for the y-Ty spectra from 100 to 857 GHz 
(presented in Fig. 12) coming from a null test (Case 1), using 
only the tSZ spectra (Case 2), and when considering both tSZ 
and tSZ-CIB spectra (Case 3). The adjustment is performed in 
the multipole range 200 > { > 2500 with 60 degrees of freedom 
(approximately 10 degrees of freedom per frequency). Thus each 
X^ value should be considered to be associated with Vdof - 10. 





V [GHz] 



100 

143 

217 

353 

545 

857 

Case 1 . . . . 

. . 600.0 

617.3 

7.9 

368.2 

172.0 

45.0 

Case 2 ... . 

4.5 

8.3 

6.8 

6.6 

41.1 

43.1 

Case 3 . . . . 

4.5 

7.8 

9.1 

8.0 

6.9 

7.0 
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5. The total tSZ-CIB cross-correlation 

In this section we investigate the all-sky tSZ-CIB cross¬ 
correlation using two different approaches: (i) the tSZ-CIB 
cross-power using a tSZ Compton-parameter map (Sect. 5.1); 
and (ii) the tSZ-CIB cross-power from a study of cross-spectra 
between Pfanck frequencies (Sect 5.2). 

5.1. Constraints on the tSZ-CIB cross-correlation from tSZ 
y-map/frequency maps cross-spectra 

This section presents the tSZ-CIB estimation using the cross¬ 
correlation between a Planck tSZ map and Planck frequency 
maps. Since the tSZ map contains CIB residuals, we carefully 
modelled these residuals in order to estimate the contribution 
from the tSZ-CIB correlations. 


_o 

(n 

E 


m 
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5.1.1. Methodology 

We compute the cross-power spectra between the Planck fre¬ 
quency maps and a reconstructed y-map"^ derived from compo¬ 
nent separation (see Planck Collaboration XXII 2015, and refer¬ 
ences therein). We choose the MILCA map and check that there 
are no significant differences with the NILC map (both maps 
are described in Planck Collaboration XXII 2015). This cross¬ 
correlation can be decomposed into four terms: 

cY^ = g{v)CY + -H g{v)CY^'^ + cf (23) 

where yciB is the CIB contamination in the tSZ map. We com¬ 
pute the uncertainties as 


Figure 13. Predicted CIB transmission in the tSZ y-map for a 
hypothetical CIB emission of 1 Kcmb at 545 GHz and at red- 
shift z. The shaded area represents the 68% confidence region, 
for uncertainties of AT^ = 2K and Afid = 0.1 in the modified 
blackbody parameters of the hypothetical CIB emission. 


where 7’‘^®(v) is the CIB emission at frequency v. Using the 
weights w,_p_v, and considering the CIB luminosity function, it 
is possible to predict the expected CIB leakage as a function of 
the redshift of the source by propagating the SED through the 
weights that are used to build the tSZ map. 


f-'y.y ^-.y.Ty ^y.T^ 

,r^Jy { { t 


(24) 


where /sky is the fraction of the sky that is unmasked. We bin 
the cross-power spectrum and deconvolve the beam and mask 
effects, then we propagate uncertainties as described in Sect. 4.2. 

The cross-correlations of a tSZ-map built from component- 
separation algorithms and Planck frequency maps are sensitive 
to both the tSZ auto-correlation and tSZ-CIB cross-correlation. 
But this cross-correlation also has a contribution produced by the 
CIB contamination to the tSZ map. In particular, this contami¬ 
nation is, by construction, highly correlated with the CIB signal 
in the frequency maps. 



5.1.2. Estimation of CIB leakage in the tSZ map 

The tSZ maps, denoted f in the following, are derived using 
component-separation methods. They are constructed through a 
linear combination of Planck frequency maps that depends on 
the angular scale and the pixel, p, as 

^M2i,p,vT/,p(v). (25) 

i.j.V 

Here Tip{v) is the Planck map at frequency v for the angular 
filter i, and are the weights of the linear combination. Then, 
the CIB contamination in the y-map is 

yciB^Y]wi,p(y)Yl^(y), (26) 

i.j.V 

* Available from http: //pla. esac. esa. int/pla/ . 


Figure 14. Expected contribution to the tSZ power spectrum for 
the true tSZ signal (black curve), for CIB leakage (red curve), 
and for tSZ-CIB leakage contribution (blue curve). The dotted 
line indicates a negative power spectrum. 


In Eig. 13, we present the expected transmission of CIB 
emission in the Planck tSZ map for a 1 Kcmb CIB source at 
545 GHz at redshift z, based on the fiducial model for the scal¬ 
ing relation presented in Sect. 2. The intensity of CIB leakage in 
the tSZ map is given by the integration of the product of CIB 
transmission (presented in Eig. 13) and the CIB scaling rela¬ 
tion at 545 GHz. Error bars account for SED variation between 
sources. In this case we assume an uncertainty of AT^ = 2 K and 
APd = 0.1 on the modified blackbody parameters. We observe 
that the CIB at low z leaks into the tSZ map with only a small 
amplitude, whereas higher-redshift CIB produces a higher, dom- 
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Figure 15. From left to right and top to bottom: observed cross-correlation between the MILCA y-map and the Planck frequency 
maps from 30 to 857 GHz. In blue are the data points, in black the CIB-cleaned data points; the red solid line is the predicted 
signal from tSZ only and the green line is the total expected signal from the tSZ signal and the tSZ-CIB correlation. All spectra are 
presented in Compton-parameter units. 


inant, level of leakage. Indeed, ILC-based component-separation 
methods tend to focus on Galactic thermal dust removal, and 
thus are less efficient at subtracting high-z CIB sources that have 
a different SED. 

The CIB power spectra have been constrained in previous 
Planck analyses (see, e.g., Planck Collaboration XVIII 2011; 
Planck Collaboration XXX 2014), as presented in Sect. 2.5. We 
can use this knowledge of the CIB power spectra to predict the 
expected CIB leakage, tcibj in the tSZ map, y. We performed 
200 Monte Carlo simulations of tSZ and CIB maps that follow 
the tSZ, CIB, and tSZ-CIB power spectra. Then, we applied to 
these simulations the weights used to build the tSZ map. Finally, 
we estimated the CIB leakage and its correlation with the tSZ ef¬ 
fect. The tSZ map signal, y, can be written asy = y + yciB- Thus, 
the spectrum of the tSZ map is + ICy™. 

In Fig. 14, we present the predicted contributions to the tSZ 
map’s power spectrum for tSZ, CIB leakage, and tSZ-CIB leak¬ 
age. We observe that at low i (below 400) the tSZ signal domi¬ 
nates CIB leakage and tSZ-CIB leakage contamination, whereas 
for higher (the signal is dominated by the CIB leakage part. The 
tSZ-CIB leakage spectrum (dotted line in the figure) is negative, 
since it is dominated by low-z (z ^ 2) CIB leakage. 

We also estimate the uncertainties on using 

the uncertainty on the CIB correlation matrix from Planck 
Collaboration XXX (2014). We derive an average uncertainty 
of 50 % on the CIB leakage amplitude in the tSZ map. This un¬ 
certainty is correlated between multipoles at a level above 90 %. 
Consequently, the uncertainty on the CIB leakage in the tSZ-map 
can be modelled as an overall amplitude factor. 

5.1.3. Results 

By cross-correlating the simulated CIB leakage signal with the 
simulated CIB at each frequency, it is also possible to predict the 
CIB leakage in the cross-spectra between tSZ map and Planck 


frequency maps. We can correct spectra (Eq. 23) using the 

estimated tSZ-CIB leakage cross-correlation term, 

and the CIB-CIB leakage cross-coiTelation term, gjy_ 

ing 

cf= cf' - g(v)C^’""“ - (27) 

Thus, the only remaining contributions are from the tSZ auto¬ 
correlation and the tSZ-CIB cross-correlation. We also propa¬ 
gate the associated uncertainties. 

Fig. 15 shows the cross-correlation of the tSZ-map and 
Planck frequency maps after correcting for the effects of the 
beam and mask, and for the terms g{v)CY'^"^ and 
As was the case for Fig. 12, all spectra are displayed in tSZ 
Compton-parameter units, and the uncertainties present a high 
degree of correlation from one frequency to another. For each 
cross-spectrum we adjust the amplitude, Atsz-ciB, of the tSZ- 
CIB contribution through a linear fit. The results of the fit are 
listed in Table 3. We obtain a maximum significance of 2.3 cr at 
857 GHz and the results are consistent with the fiducial model. 

Table 3. Best-fit values for the tSZ-CIB amplitude, Atsz-ciB, 
using the fiducial model as reference. 


V [GHz] 

Atsz-CIB 

AAtsz-ciB 

100. 

-3.6 

3.8 

143. 

-1.6 

3.7 

353 . 

2.0 

2.0 

545 . 

1.7 

1.4 

857 . 

1.6 

0.7 


The combined constraints from 353, 545, and 857 GHz mea¬ 
surements, as well as the covariance structure of the measure¬ 
ment uncertainties, yield an estimate of Atsz-ciB = 1.3 + 0.4. 
The uncertainties are dominated by CIB leakage subtraction. 
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which leads to highly correlated uncertainties of the 353, 545, 
and 857 GHz spectra. The optimal linear estimator we use to 
calculate Atsz-ciB probes the measurements for the known fre¬ 
quency trend of the CIB leakage in order to correct for this. 
Since the CIB leakage has affected all three of these correlated 
measurements in the same direction, the combined constraint of 
1.3 + 0.4 is below the range of the individual (uncorrected) con¬ 
straints ranging from 1.6 to 2.0. Table 3 contains the necessary 
covariance information used here. 

5.2. Constraints on tSZ-CIB cross-correlation from Planck 
frequency maps 

As a last approach, we explore the direct cross-correlation be¬ 
tween Planck frequency maps. 

5.2.1. Methodology 

In terms of tSZ and CIB components the cross-spectra between 
frequencies v and V can be written as 

cf ^g{v)g{V)Cy + 

-I- H- C°*“(v), (28) 

where Cf*“(v) accounts for the contribution of all components 
except for tSZ and CIB. We compute the cross-spectra between 
Planck frequency maps from 100 to 857 GHz as 

cy'^+cy'' 


where subscripts 1 and 2 label the “half-ring” Planck maps. This 
process allows us to produce power spectra without the noise 
contribution. We also compute the covariance between spectra 
as 


cov(c;’",cr) 



'ey + 




4(2f’H- 

l)/sky 



Cf""' + 




4(2^ H- 

l)/sky 



cy + 


'r^vy 


4(2^ H- 

l)/sky 




^V2,y;" 



't ^ 

4(2e + l)/,ky 


(30) 


We correct the cross-spectra for beam and mask effects, using 
the same Galactic mask as in Sect. 4.2, removing 60% of the 
sky, and we propagate uncertainties on cross-power spectra as 
described in Sect. 4.2. 

The tSZ and CIB contributions to C^’are contaminated 
by other astrophysical emission. We remove the CMB con¬ 
tribution in using the Planck best-ht cosmology (Planck 
Collaboration XIII 2015). We note that the Planck CMB maps 
suffer from tSZ and CIB residuals, so they cannot be used for 
our purpose. 


5.2.2. Estimation of tSZ-CIB amplitude 

We ht thermal dust, radio sources, tSZ, CIB (that accounts for 
the total fluctuations in extragalactic infrared emission), and 


tSZ-CIB amplitudes, Adu.st, Arad, Atsz, Acib, and Atsz-ciB, re¬ 
spectively, through a linear ht. For the dust spectrum we as¬ 
sume Cf oc ^ 2.8 (Planck Collaboration XXX 2014), and for 
radio sources Ci oc For the tSZ-CIB correlation, tSZ, and 
CIB power spectra, we use templates computed as presented in 
Sect. 2.3. This gives us 

CY = Atsz.g(v).g(v')C^’^ 

+ Arsz-CIB [g(y)cf +g(v')C^’°"'’''] 

+ Adust/dust(v)/dust(v')^~^*^ 

+ Arad/rad(v)/rad(v')- (31) 

Here /dust and /ad give the frequency dependence of thermal 
dust and radio point sources, respectively. For thermal dust we 
assume a modihed blackbody emission law, with /d = 155 
and Td = 20.8 K (Planck Collaboration XI 2014). For radio 
point sources we assume a spectral index = -0.7 (Planck 
Collaboration XXVIII 2014). The adjustment of the amplitudes 
and the estimation of the amplitude covariance matrix are per¬ 
formed simultaneously on the six auto-spectra and 15 cross¬ 
spectra from ^ = 50 to f = 2000. 

For tSZ and CIB spectra we reconstruct amplitudes compat¬ 
ible with previous constraints (see Sect. 2.5): Acib = 0.98 + 0.03 
for the CIB; and Atsz = 1-01 + 0.05 for tSZ. For the tSZ-CIB 
contribution we obtain Atsz-ciB = 1.19 + 0.30. Thus, we obtain 
a detection of the tSZ-CIB cross-correlation at 4cr, consistent 
with the model. In Fig. 16, we present the correlation matrix 
between cross-spectra component amplitudes. The highest de¬ 
generacy occurs between tSZ and tSZ-CIB amplitudes, with a 
correlation of -50%. We also note that CIB and radio contri¬ 
butions are signihcantly degenerate, with tSZ-CIB correlation 
amplitudes of -28 % and 29 %, respectively. 



0.5 


- 1.0 


^ISZ-ClB -^ClB -^tSZ -^rad -^dust 

Cross-spectra components 


Figure 16. Correlation matrix for multi-frequency cross-spectra 
components from Eq. (31). 
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6. Conclusions and discussion 

We have performed a comprehensive analysis of the infrared 
emission from galaxy clusters. We have proposed a model of 
the tSZ-CIB correlation derived from coherent modelling of 
both the tSZ and CIB at galaxy clusters. We have shown that 
the models of the tSZ and CIB power spectra reproduce fairly 
well the observed power spectra from the Planck data. Using 
this approach, we have been able to predict the expected tSZ- 
CIB cross-spectrum. Our predictions are consistent with previ¬ 
ous work reported in the literature (Addison et al. 2012; Zahn 
et al. 2012 ). 

We have demonstrated that the CIB scaling relation from 
Planck Collaboration XXX (2014) is able to reproduce the ob¬ 
served stacked SED of Planck confirmed clusters. We have also 
set constraints on the profile of the this emission and found that 
the infrared emission is more extended than the tSZ profile. We 
also find that the infrared profile is more extended than seen in 
previous work (see e.g., Xia et al. 2012) based on numerical sim¬ 
ulation (Bullock et al. 2001). Fitting for the concentration of an 
NEW profile, the infrared emission shows C 500 = 1-OO^lio [y This 
demonstrates that the infrared brightness of cluster-core galaxies 
is lower than that of outlying galaxies. 

We have presented three distinct approaches for constrain¬ 
ing the tSZ-CIB cross-correlation level: (i) using confirmed tSZ 
clusters; (ii) through cross-correlating a tSZ Compton parameter 
map with Planck frequency maps; and (iii) by directly cross- 
correlating Planck frequency maps. We have compared these 
analyses with the predictions from the model and derived consis¬ 
tent results. We obtain: (i) a detection of the tSZ-IR correlation 
at 6 cr; (ii) an amplitude Atsz-ciB = 1-5 + 0.5; and (iii) an ampli¬ 
tude Atsz-ciB = 1.2 + 0.3. At 143 GHz these values correspond 
to correlation coefficients at f = 3000 of: (ii) p - 0.18 + 0.07; 
and (iii) p - 0.16 + 0.04. These results are consistent with pre¬ 
vious analyses by the ACT collaboration, which set upper limits 
p < 0.2 (Dunkley et al. 2013) and by the SPT collaboration, 
which found p - O.IUqq® (George et al. 2014). 

Our results, with a detection of the full tSZ-CIB cross¬ 
correlation amplitude at 4 cr, provide the tightest constraint so far 
on the tSZ-CIB correlation factor. Such constraints on the tSZ- 
CIB cross-correlation are needed to perform an accurate mea¬ 
surement of the tSZ power spectrum. Beyond power spectrum 
analyses, the tSZ-CIB cross-correlation is also a major issue for 
relativistic tSZ studies, since CIB emission toward galaxy clus¬ 
ters mimics the relativistic tSZ correction and thus could produce 
significant bias if not accounted for properly. This 4 cr measure¬ 
ment of the amplitude of the tSZ-CIB correlation will also be 
important for estimates of the “kinetic” SZ power spectrum. 
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